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Abstract 


In  this  paper  we  study  Che  divergence  operator  acting  on  continuous 


piecewise  polynomials  of  degree  p  +  1  ,  pi  3  ,  on  triangulations  of  a  plane 

J  _ 

polygonal  domain  JSC  .  -WS  give  a  characterization  of  the  range  of  the  diver¬ 


gence  operator  and  the  full  details  of  a  combinatorial  verification  of  this.  t 
As  the  central  result  we  show  that  for  very  general  families  of 


meshes  it  is  possible  to  find  a  maximal  right  inverse  for  the  divergence 
operator  with  a  SCL^;  H^)  norm  which  is  bounded  independently  of  the  mesh 
size.  The  norm  of  this  right  inverse  grows  at  most  algebraically  with  p  , 
but  it  necessarily  blows  up  as  a  certain  measure  of  singularity  of  the  meshes 


approaches  0  . 


1.  Introduction 


Incompressibility  constraints,  such  as  constraints  on  the  divergence  of 
a  velocity  field  or  a  displacement  field,  occur  in  many  equations  of  physical 
interest,  e.g.  the  Navier-Stokes  equations  or  the  equations  of  elasticity. 

When  analyzing  the  stability  of  finite  element  approximations  to  these 
equations  a  central  question  concerns  the  behaviour  of  the  divergence  operator, 
or  a  discrete  version  thereof,  on  the  corresponding  spaces  of  piecewise 
polynomials  (see  for  instance  [4,7,17,19]).  It  is  well  documented  that  continuous 
piecewise  polynomials  of  low  degree  applied  directly  to  the  velocity-  (or 
displacement-)  formulation  are  often  inadequate,  due  to  the  lack  of  a  uniformly 
bounded  right  inverse  for  the  divergence  operator.  This  has  led  various 
authors  to  study  non-conforming  low  order  elements  in  connection  with  mixed 
formulations.  The  analysis  in  this  paper  points  in  another  direction:  our 
results  imply  that  continous  piecewise  polynomials  of  degree  four  or  higher 
directly  applied  to  the  velocity  -  (or  displacement-)  formulation  lead  to 
optimal  (uniform)  convergence  rates  (for  a  discussion  of  this,  see  [12]). 

The  paper  [18]  contains  a  characterization  of  the  range  of  the  divergence 
operator  on  spaces  of  continuous  piecewise  polynomials  of  degree  p  +  1  ,  p  *  3  , 
on  an  arbitrary  triangulation  (Theorem  2.1  and  Remark  2.1);  it  also  gives  a 
proof  of  the  fact  that  on  a  fixed  triangulation  it  is  always  possible  to 
construct  a  maximal  right  inverse  for  the  divergence  operator,  the  norm  of 
which  grows  at  most  algebraically  with  p  (Theorem  2.1).  These  results  were 
used  to  prove  that  the  so-called  p-version  of  the  finite  element  method,  when 
applied  directly  to  the  displacement  formulation  of  plane  strain  elasticity, 
converges  at  optimal  rate  independently  of  the  value  of  Poisson’s  ratio  ([19]). 

The  analysis  presented  in  this  paper  extends  the  results  of  [18]  in 
a  rather  surprising  way  -  it  shows  that  the  aforementioned  right  inverse 
(p  >,  3)  has  a  Bd^jH*)  operator  norm  which  is  bounded  Independently  of  the 


2 


mesh  size  of  the  triangulation.  This  uniform  bound  can  only  hold  provided  a 
certain  measure  of  singularity  of  the  meshes  is  bounded  away  from  zero  (cf. 

Example  3.1).  Available  numerical  experiments  (cf.  [16])  and  recent  theoretical 
results  (cf.  [12])  indicate  that  a  similar  bound  does  not  exists  for  p  <  3  . 

For  reasons  of  exposition  we  have  chosen  to  express  our  main  result  in 
terms  of  a  bound  for  the  norm  of  a  maximal  right  inverse  for  the  divergence 
operator.  It  is  easy  to  see  (cf.  section  5)  that  this  is  equivalent  to  a 
uniform,  positive  lower  bound  for  the  expression 

7  L  T-2  *  d-/\l\\ .1  11*11, 

—  H  L 

as  studied  by  other  authors  (here  V  varies  over  the  space  of  piecewise 
polynomials  and  <j>  varies  over  the  divergence  of  this  space). 

The  organization  of  this  paper  is  as  follows:  in  section  2  we  introduce 
the  necessary  notation  concerning  the  triangulations  and  the  polynomial 
subspaces.  It  should  be  emphasized  that  our  triangulations  are  quite  general 
and  only  restricted  by  the  assumption  of  quasiuniformity.  Section  3  independ¬ 
ently  characterizes  the  range  of  the  divergence  operator  acting  on  continuous 
piecewise  polynomials  of  degree  p  +  1  ,  p  >,  3  .  Combinatorial  proofs  are  carried 
out  both  with  and  without  boundary  conditions;  in  the  latter  case  the  argu¬ 
ment  is  identical  to  one  found  in  [18]  and  depends  crucially  on  the  formula 
for  the  dimension  of  piecewise  polynomials  proven  in  [10];  in  the  first 

case  we  have  to  establish  a  similar  formula  for  piecewise  polynomials 

that  vanish  to  second  order  on  the  boundary  (this  is  done  in  section  6). 

Section  4  and  5  contain  the  proof  of  the  main  theorem,  the  existence  of  a 
uniformly  bounded  maximal  right  inverse.  The  analysis  relies  heavily  on  [18], 
but  an  important  new  element  is  the  localization  procedure  formulated  in 
lemmas  4.3,  5.1  and  5.2.  The  idea  behind  lemmas  5.1  and  5.2  is  in  many 


ways  similar  to  that  underlying  the  macro-element  technique  and  the  correspond¬ 
ing  local  test  for  stability  found  in  [3]  or  [14].  Much  of  the  rest  of  the 
proof  of  the  main  theorem  consists  of  verifying  that  the  constants  in  various 
of  the  estimates  found  in  [18]  scale  appropriately  with  the  mesh  size. 

The  attention  in  this  paper  is  restricted  to  plane  domains;  it  should 

3 

be  interesting  to  see  if  a  similar  analysis  could  be  carried  out  in  E 
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2 .  Notation 

2 

Throughout  this  papor  donor  os  a  hounded  polvgon.il  domain  in  K 

£  =  ^jTijN(h)  ^  0  <  h  .<  1  ,  is  a  family  of  triangulations  of  ft  ,  para- 

h  1  i”1  _h 

metrized  by  mesh  size  h  .  To  be  more  precise:  the  i  ^  ,  1  {  i  $  N(h)  , 

for  fixed  h  ,  are  disjoint  triangles  with 

diam  $  h 


and 


N(h)  — r 

U 

i=l 


SI  . 


An  edge  of  a  triangle  of  ^  is  called  an  internal  edge  of  £  if  its  interior 
lies  in  ft  (not  on  Ift  )  .  We  assume  that  no  vertex  of  a  triangle  of 
falLs  in  the  interior  of  an  internal  edge  of  £  .  This  does  not  prevent  boundary 
edges  from  having  vertices  in  their  interior  (as  in  Fig.  2  and  Fig.  4). 

Furthermore  we  assume  that  the  family  ,  0  <  h  '  1  ,  is  quasiuniform  in 

I  III’  sen.'. I '  l  ll.it 

(2.1)  0Qh  <  0(7)  vr  €  lh  ,  0  <  h  i  l  , 


where  p(T)  denotes  the  supremum  of  diameters  of  discs  contained  in  T  ,  and 
0  ■  Pq  •  In  the  rest  of  this  section  and  all  of  the  next  we  shall,  to  simplify 
notation,  omit  the  subscript  h  when  referring  to  a  fixed  triangulation. 

If  is  an  arbitrary  subset  of  the  triangulation  £  ,  we  then  define 

the  corresponding  polygonal  domain* 

(2.2)  ft(I')  *  interior  (  U  THft)  . 

For  any  integer  p  >,  0  and  r  =  0  or  1 


p(p]’r(£') 


*0(1  in  the  definition  of  ft( £  )  matters  only  when  interior  (ft)  4  ft  ,  e.g.  when 
ft  is  a  slit  doma in . 


& 


'  ( 
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denotes  the  set  of  functions  in  Cr(Sl(£  ))  that  are  given  by  a  polynomial 

n  f 

of  degree  <  p  on  each  of  the  triangles  of  2 

An  internal  vertex  of  [  is  a  vertex  that  lies  in  Q([  )  (not  on 
3ft(£'))  .  We  shall  say  that  an  internal  vertex  of  [  is  singular  if  the 
edges  meeting  at  this  vertex  fall  on  two  straight  lines  (cf.  [10]). 


Fig.  1.  Singular  internal  vertex  x^ 


Following  [18]  we  introduce,  for  p  >  0  ,  the  space 


of  functions  ,  i  ,  which  are  given  by  a  polynomial  of  degree  $  p  on  each 

individual  triangle  (no  continuity  requirements)  and  which  have  the  property  that 

-* 

R.l :  at  any  singular  internal  vertex  of  2  *  *n  * 


1  i 

2  (xq)  ■  0 

i-i  1  ^ 


where  ^(x^)  ■  (Xq)  and  7^...,7^  are  the  triangles  meeting 
at  Xq  *  numbered  consecutively,  as  shown  in  Fig.  1. 

At.  explanation  for  the  requirement  (R.l)  is  most  easily  given  by  the 
following  simple  observation. 


,1 


6 


Proposition  2 . 1 

For  any  £  £  and  any  p  i  0  , 

7.(p[p+l],°(£*)  x  p(p+l]»0(£»)}  £  pIP1*"1^')  . 

The  proof  of  this  proposition  consists  of  a  straightforward  calculation, 
the  details  of  which  we  le^ve  as  an  exercise.  Special  cases  of  this  result 
have  been  used  by  other  authors,  e. g.  Mercier  [9]  and  Fix  et  al.  [6]. 

When  homogeneous  Dirichlet  boundary  conditions  are  imposed,  a  new  set  of 
requirements  become  important.  Let 

Ptpl,r(£’>  *  P  *  r  +  1  ,  r  -  0,1 

denote  the  subspace  of  P^,r(£')  consisting  of  those  functions  that  vanish 

O  r  1 

to  order  r  +  1  on  3fl(£')  ;  that  is,  functions  in  P  P  ,r(£')  are  always 
zero  on  3Q([')  ,  and  in  addition,  functions  in  P^’1^')  are  required  to 
have  a  vanishing  normal  derivative. 

Remark  2 . 1 

In  this  paper  we  use  the  very  natural  convention,  that  a  point  on  an(T')  , 
which  is  a  vertex  for  k  different  parts  of  30(2')  .  be  considered  k  different 
boundary  vertices.  As  an  example  there  are  two  different  boundary  vertices  at  the 
point  P  in  Fig.  2.  A  similar  convention  is  applied  to  edges  that  lie  on  "internal" 
boundaries.  These  are  considered  two  different  boundary  edges  if  they  are 
common  to  two  different  triangles  of  .  Note  that,  conforming  with  this 

convention,  our  definitions  of  piecewise  polynomial  spaces  do  not  impose  any 
continuity  conditions  at  vertices  or  edges  where  the  boundary  intersects  itself. 

□ 


4.  " 


f 


Point  at  t 

Afferent  ls  “■»»«.«  to  be  too 

A  vertex  3{1(T’)  t 

311  ^e  edges  of  T  d  3  S±n8Ular  bo^dary  vertex  of  V  if 

There  are  four  ' =«•  vertex  «,  _  ^  /  " 

*>«  Possible  ==ofl8uratw  “™i8ht  Uoea. 

^  3.  (The  toertb  caw  ^  S1"8U1"  vertex,  aa 

(l81  •*»«  it  ais„  Uluatr.tea  the  .  ''  '  d"f"8  «-t  to 

ln  the  lnt"l°r  <*  •  boundary  edge  T  ’  °'  *  b°'"daCi'  V"'«  ^ 
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We  let 


pE^-V)  ,  p  *  o 

denote  the  subspace  of  P^’  )  consisting  of  functions,  $  ,  which 

additionally  satisfy  the  following  two  requirements,  that 

r » 

R. 2 :  At  any  singular  boundary  vertex  of  l  ,  , 

k 

l  <-i)Yoo  =  o  , 

i=i  1  ^ 

where  ^(Xq)  =  4>  I  y  »  and  , . . . , are  the  triangles  of 

r-  '  ^ 

l  meeting  at  j£q  (k  can  be  any  number  from  1  to  4,  and  the 
triangles  are  numbered  consecutively  as  shown  in  Fig.  3). 

R.  3 :  For  any  connected  component  of  £2(J  )  ,  £2"  , 

$  dx  =  0  . 

It  is  a  simple  exercise  to  show  that  the  following  holds. 

Proposition  2 . 2 

For  any  £  c  £  and  any  p  >  0  , 

v.(p[P+l],0(^)  x  p[ p+1  ] , 0 (£ ' } }  £  p[p],-l(^}  # 

Our  notation  for  Sobolev  spaces  is  standard:  if  £2 '  £  £2  is  a  polygonal 
(sub)  domain  then  H  (£2')  ,  k  €  IN  ,  denotes  the  set  of  functions  with  deriva¬ 
tives  of  order  S  k  in  L0  (£2  ’ )  ;  the  corresponding  norm  is  denoted  II 'II.  0i 
H^(£2')  is  the  closure  of  C”(£2')  in  H^(i2')  . 


For  Che  analysis  of  finite  element  discretizations  of  equations  with  a 
divergence  constraint  it  is  important  to  have  precise  information  about  the 
range  of  the  divergence  operator  on  the  finite  dimensional  subspaces.  In 
general  a  uniform  norm  estimate  of  a  right  inverse  is  sufficient  to  guarantee 
stability,  however,  in  order  to  estimate  the  convergence  rate,  the  algebraic 
character  of  the  range  of  the  divergence  operator  has  to  be  known.  In  the 
present  situation  it  furthermore  turns  out  that  the  characterization  of  the 
range  automatically  leads  to  a  necessary  condition  for  the  existence  of  a 
uniformly  bounded  right  inverse.  For  ease  of  notation  we  omit  the  subscript 
h  when  referring  to  a  fixed  triangulation.  The  following  result  was 
proven  in  [18]. 

Proposition  3.1 

For  any  [  £  [  and  any  p  >,  3  the  divergence  operator  maps 

p(p+l],0(2')  x  p[p+l].0(J') 

onto 

plpl. _ 

As  shown  in  [18]  this  result  permits  a  simple  combinatorial  proof.  We  give 
the  full  details  of  the  combinatorial  argument  below. 

Consider  first  the  case  that  R(£  )  is  simply  connected.  The  curl  operator 

7  ‘  *  •  (k;  ♦  •  -  %  ♦) 

[p+2],l 


maps  P 


(£  )  onto  the  nullspace  of  the  divergence  operator 
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NP+V.)  c  PtP+1]*°(^)  X  P[ P+1  . 

In  [10]  It  is  shown,  that  for  p  >,  3 

(3.1)  dim  (P[p+2],1(£'))  = 
l/2(p+3)(p+4)T  -  (2p+5)EQ  +  3VQ  +  aQ  , 

where  T  is  the  number  of  triangles,  Eg  is  the  number  of  internal  edges, 
Vq  is  the  number  of  internal  vertices  and  Og  is  the  number  of  singular 
internal  vertices,  all  of  the  subset  £  £.  I  •  Since  the  nullspace  of  the 

curl  operator  consists  of  only  the  constants,  it  follows  from  Grassmann's 
dimension  formula,  that 

(3.2)  dim(Np+1(7-))  =  dim(P[p+2J,1(£’))  -  1  . 

If  RP ( V • )  denotes  the  range  of  the  divergence  operator  acting  on 
p[p+l]i0(£')  x  p[p+l].0(J')  ^  then  the  same  dimension  formula  gives  that 

(3.3)  dim(RP(V-))  =  dim(P[p+1]’°(f )  *  Ptp+1 1 ’°(^ ' ) ) 

-  dim(NP+1 (V • ) )  . 

The  first  term  in  the  right  hand  side  of  (3.3)  is  easily  found  to  be 

(3.4)  (p-l)pT  +  2pE  +  2V  , 

where  T  is  as  before,  E  denotes  the  total  number  of  edges  and  V  the 
total  number  of  vertices  of  £  .  Inserting  (3.1),  (3.2)  and  (3.4)  into 


(3.3)  one  gets 


(3.5) 


dira(RF (V • ) )  =  1/2 (p  -9p-12 )T  +  2(p+l)(E+E0> 


-  2E  +  3E0  +  2V  -  3Vq 


ao 


+  1 


=  1/2 (p+3)pT  +  E  -  V  -  oQ  +  1  , 


with  the  second  identity  based  entirely  on  the  relations  V  -  *  E  -  E^  and 

E  +  Eq  =  3T  .  Euler's  formula  states  that 

T  -  E  +  V  =  1  , 

and  in  combination  with  (3.5)  this  gives  1 


(3.6)  dim(RP(V • ) )  =  1/2 (p+2) (p+l)T  -  oQ  . 

The  right  hand  side  of  (3.6)  is  exactly  the  expression  for  the  dimension  of 
p[pl>  1(£')  ,  This  observation  together  with  Proposition  2.1  implies  that 

Rp ( V • )  =  pip] ,-i(^»)  . 

if  n(I')  is  not  simply  connected  then  we  extend  any  function  in  P^p^’  ^ 
by  piecewise  linear  functions  onto  triangles  filling  the  holes  of  ft(£')  •  Thi 
can  be  done  in  such  a  way  that  the  extension  is  still  in  P^p^’  ^  ancj  we  may 
now  rely  on  the  previous  argument  to  ensure  the  existence  of  a  field  in 
ptp+lI.O  x  p[p+l].0  w£Ch  the  right  divergence.  (It  has  here  implicitly  been 
assumed  that  the  holesof  fl(£')  have  boundaries  that  are  not  self intersecting; 
selfintersecting  boundaries  can  be  dealt  with  by  a  perturbation  argument.) 

With  homogeneous  Dirichlet  boundary  conditions  the  corresponding  result  is 

Proposition  3.2 

For  any  £  £  £  and  any  p  5  3  the  divergence  operator  maps 


p[p+l],0.r*.  n[p+l],0.r» 


onto 
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plpJ’-L<r> . 


The  analysis  given  in  [18]  verifies  this  for  p  sufficiently  large  by  an 
approximation  argument.  At  the  end  of  section  4  we  show  how  the  result  Is  obtained 
for  general  [ '  l  and  p;  3  .  Ior  completeness  we  briefly  outline  the  key 
ingredients  of  a  combinatorial  argument.  Suppose  that  )  is  simply 

connected;  in  that  case 


(3.7)  dim(P[p+21,1(r>)  = 

l/2p(p-5)T  +  (2p-l)E0  +  3VQ  +  o  , 

where  o  denotes  the  total  number  of  singular  vertices  of  £  .  The  formula 

(3.7)  is  verified  in  section  6  by  a  method  based  on  [10],  We  furthermore 
know  that 


(3.8)  dim(P[p+11’°(f )  x  Ptp+11’°(r»  “  (P-DpT  +  2pE0  +  2V0 


dim(P[p,’"1(r»  =  l/2(p+2)(p+l)T  -  a  -  1  . 


The  formulae  (3.7)  and  (3.8)  in  combination  with  an  argument  like  the  preceding 
may  now  he  applied  to  prove  Proposition  J.2  whenever  U ( ^ ' )  is  simply  connected. 
Non-simplv  connected  domains  may  he  treated  by  a  slight  variation  of  this 
argument  (cf.  Remark  b.l). 

Remark  3 . 1 

Propositions  3.1  and  3.2  remain  valid  also  for  p  =  2,1  or  0  on 
any  £'  such  that  Si(£')  is  simply  connected  and  the  formula  (3.1), 
respectively  (3.7),  holds.  The  formula  (3.1),  which  was  conjectured  by 
St  rang  [ !  r>  | ,  has  been  verified  for  certain  triangulations  J'  (in  decreasing 
generality,  as  p  decreases)  by  Morgan  and  Scott  [11].  The  formula  (3.7), 
however,  fails  on  tlu*  most  natural  triangulutions  as  soon  as  p  <  2  .  For 
a  more  detailed  discussion  we  refer  to  [12]. 
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The  main  goal  in  this  paper  is  to  verify  the  existence  of  a  maximal 
right  inverse  for  the  divergence  operator,  the  norm  of  which  is  bounded 
uniformly  in  the  mesh  size  and  grows  at  most  algebrically  with  p  .  It  turns 
out  that  our  proof  of  this  fact  does  not  depend  on  Propositions  3.1  and 
3.2,  to  the  contrary,  it  provides  independent  proofs  of  these.  However,  these 
propositions  demonstrate  the  necessity  of  a  certain  non-degeneracy  condition 
on  the  meshes,  if  one  wants  to  obtain  a  uniformly  bounded  right  inverse. 

r»  t 

Let  X{)  denote  any  non-singular  vertex  of  £  and  ^et  9^  ,  1  {  i  {  k  , 
be  the  angles  of  the  triangles  ^  ,  1  j  i  j  k  ,  meeting  at  Xq  (the  triangles 
are  numbered  consecutively  as  before).  We  define 

R(Xq)  *  max{ j 0^+0^-ff |  ;  l<i,j<k  and  i-j  »  1  mod  k)  ; 

R(Xq)  thus  measures  how  close  x^  is  to  being  singular.  We  furthermore  set 

(3.9)  R([')  =  min{R(x^)  :  x^  is  a  non-singular 

internal  vertex  of  £  } 

and 

o  T 

(3.10)  R(^  )  =  rain{R(xQ)  :  Xq  is  a  non-singular 

vertex  of  £  }  . 


Example  3. 1 

r*  (5  P  (5 

Let  i  be  the  simple  traingulation  shown  in  Fig.  4,  with  R (£  )  ■  5  . 
Let  ,  5  small,  be  the  piecewise  constant  that  is  given  by 


fl  in  T 


0  otherwise . 


/  t  / 

/ 

/  '' 

/  ' 

(' ' 

V 


Fig.  4.  I1 

[4]  0  r6 

Proposition  3.1  ensures  that  for  6  >  0  there  exists  €  P  *  Q  ) 
x  P[Al,0(£5)  With 


but  ||  V  ^  ^  ,5  cannot  stay  bounded  as  5  -*-  0  .  If  ll  V  g  |[  ^  flj  ?  C  ,  uniformly 

as  .5  0  ,  then  we  could  extract  a  weakly  convergent  subsequence,  which  would 


converge  to  a  field 


V-  «  P(a1,0(£°)  x  ptA]’°(£°)  , 


satisfying 


1  in  T, 


D  otherwise. 


rO 

This  is  a  contradiction,  since  x^  is  a  singular  internal  vertex  for  2  •  [ — | 

The  previous  example  shows  that  it  is  in  general  necessary  to  have  (3.9) 
(or  (3.10>)  bounded  from  below  In  order  to  establish  uniform  bounds  for  the 
o  ( ! . , ;  H 1 )  norm  of  im.ixim.il  right  Inverse  for  the  divergence  operator. 


i 


/ 
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4.  Local  construct  ion  of  <3  r  ight  inverse  for  V*  . 

The  first  in  a  series  of  lemmas  is  an  extension  of  Lemma  2.6  in  [18], 


Lemma  4 . 1 


Assume  that 


R(L)  5  6  >  0  , 


where  R(£^)  is  the  measure  of  singularity  introduced  in  (3.9),  and  6  is 
independent  of  h  .  Let  ^  denote  any  subset  of  ,  and  let  4>  be  any 

element  of  P^*  )  .  There  exists  V  €  )  x  such  that 


(4.1a) 


(4.1b) 


4>  -  7‘V  =  0  at  all  vertices  of  ,  and 


'-11  i  Ml'h) 


\  $  C(p+1) 


o,n(D 


with  constant  C  and  K  that  are 
independent  of  ,  h,p  and  $  . 


Proof : 


Let  and  be  two  adjacent  unit  sized  triangles,  as 


shown  in 


Fig.  5. 


Fig.  5.  Two  adjacent  triangles. 
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If  a  is  any  constant,  then  it  is  possible  to  find  a  continuous  piecewise 
cubic  field  U  on  T^UT^  ,  satisfying 

V»£  =  a  at  Xq 

(4.2)  V*_U  =  0  at  all  other  vertices,  and 

U  =  0  on  scr^uT^)  . 

From  the  construction  in  [18]  it  follows  that 

!|u||  i  t  iff  ^  ^ I a I 

A*  lu  2 

where  C  only  depends  on  the  minimal  angle  of  T  ,  i  *  1,2  .  If  furthermore 
0^  +  0^  T6  it  and  a^  ,  a^  are  any  two  constants,  then  one  can  find  a  continuous 
piecewise  cubic  field  U'  on  T^UT ^  ,  such  that 

v‘E'|t  (Xq)  =  a j:  for  1  =  l.2. 

(4.3)  V*u'  =0  at  all  other  vertices 

U  *  0  on  a  (T  UT  ) 


where  C  depends  on  the  minimal  angle  of  T  ,  i  =  1,2  ,  and  j |  . 

Let  Xq  be  a  non-singular  internal  vertex  with  N  corresponding  triangles 
of  unit  size,  and  let  a^  ,  1  $  i  $  N  ,  be  N  arbitrary  constants.  Using 
(4.2),  (4.3)  and  the  same  argument  as  in  [18]  we  obtain  a  continuous  piecewise 
cubic  field  W  with 


(4.4) 


v’Wjy  (xq)  =  a^  for  1  <  i  $  N  , 

7»W  “0  at  all  other  vertices 


W  *  0  on 
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This  field  can  be  estimated  by 

N 

(4*5)  II  W||  nr  <  C  £  I  a .  |  , 

’  1  i=l 

where  C  only  depends  on  the  minimal  angle  of  T  ,  1  $  i  (  N  ,  and  R(Xq) 

(C  blows  up  when  either  of  these  become  too  small).  At  any  singular  internal 

vertex  we  may  similarly  find  a  continuous  piecewise  cubic  field  satisfying 

4 

(4.4),  (4.5)  provided  £  (-1)  a.  «  0  .  The  constant  C  here  depends  only  on 

i=l  1 

the  minimal  angle.  Since  we  are  not  imposing  any  boundary  conditions  (4.4) 
and  (4.5)  can  also  be  satisfied  for  any  boundary  vertex  and  any  set  of 
constants  a^  ,  with  a  constant  C  that  only  depends  on  the  minimal  angle. 

By  rescaling  we  see  that  all  these  versions  of  (4.4),  (4.5)  remain  valid 
with  a  constant  that  is  Ch  ,  where  h  is  the  size  of  the  triangles.  For 
each  vertex  Xq  Ih  we  select  a.  ,  1  .<  i  s  N  ,  to  be  <f> 

previous  construction  then  leads  to 

||w||  <  c(p+i)K|M| 

i.'jTv  o.t/r 

i  i 

for  K  >  2  (of.  [18  |).  Adding  the  individual  W’s  we  arrive  at  a 

* 

field  V  ,  satisfying  (4.1a)  and  (4.1b)  .  The  constant  C  is  independent  of 
since  both  R(£^)  anc*  minimal  angle  are  bounded  away  from  0  (the 
latter  because  of  the  quasiuniformity  assumption).  | _ ) 

Remark  4 . 1 

Assume  that  R(^)  5  6  >  0  and  that  $  €  P^’  *■(£  )  with  «  0  at 
the  boundary  vertices  of  fl(£  )  .  Then  it  is  possible  to  find  V  €  P^’^(£^) 
x  )  such  that  (4.1a-b)  hold.  It  is  crucial  that  $  ■  0  at  the 

vertices  on  3fi(£^)  provided  we  want  to  maintain  R(£^)  39  c^e  measure  of 

O  f 

singularity.  If  we  make  the  alternate  assumption  that  R(£^)  5  S  >  0  then 


|Th(^o)  ;  the 
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it  is  possible  to  find  V  €  )  x  satisfying  (4.1a-b)  for 

any  f>  €  P^’  .  These  slight  variations  of  Lemma  4.1  follow  by  a  proof 

very  similar  to  the  previous.  |— | 

Let  be  two  arbitrary  triangles  of  ^  ,  with  a  common  edge 

(as  in  Fig.  6). 


\ 


/ 


\  / 
\  / 

\  / 


/ 


Fig.  6. 


Denote  by 


*J(«) 


h  ,  „h 

a.x  +  8 . x  +  y  =  0 
i  1  i  2  i 


1  S  i  .<  4  , 


the  four  lines  on  which  the  remaining  edges  lie,  and  define 


£2£2 

12 


x  €  T 


h 

1 


ip(x) 


x  €  r: 


where  c  is  chosen  such  that  i>  is  continuous  in  UT^  .  Let  n  be  a 

normal  direction  to  the  common  edge,  and  introduce 
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WOO  =  d  i|i(x)  ji  . 


Any  such  W  satisfies 


V*W  dx  =  -  |  7«W  dx 


and  by  choosing  d  ^  0  appropriately  we  thus  obtain 


Lemma  4 . 2 

rh  ,  -rh 


Let  and  be  two  triangles  of  with  a  common  edge.  It  is 

possible  to  find  a  continuous  field  W  such  that 


(4.6a)  W  is  given  by  polynomials  of  degree  $  4  on  each  of  the 
triangles  ,  and  W  =  0  on  SCT^UT^)  , 


(4.6b)  7*W  =0  at  all  vertices  of  7\  ,  i  =  1,2, 


(4.6c) 


-[ 


V*W  dx  = 


V-W  dx  =  1  , 


(4.6d)  1 1 W 1 1  _  .  <  Dh  ^  where  D  is  independent  of  and  h  . 


Note 


:  In  the  estimate  (4.6d)  we  have  used  the  fact  that  the  triangulation 


satisfies  a  minimal  angle  condition  due  to  the  assumption  of  quasiuniformity. 


Def init ion 


f  h  / 

A  subset  >,  =  {Tnj  of  7,  is  called  connected  if  the  corresponding 

Lh  l  Lh  „ 

r*  *  — 

polygonal  domain  SKK  )  =  interior(  U  T  flfi)  is  connected. 

"  i-1  1 


Lemma  4 . 3 

Let  D  be  the  same  constant  as  in  the  previous  lemma.  For  any  connected 

i  h  £  £ 

subset  7^  *  £  7^  and  any  set  of  numbers  ,  W*' 
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l 


b.  =  0  , 
1  * 


one  can  find  \/  €  )  x  pIM»0^£  j  satisfying 

(4.7a)  V*V  =  0  at  all  vertices  of  ^  , 


(4.7b) 


V*V  dx  =  b  ,  1  j:  i  $  £  ,  and 


(4.7c) 


W  1  ,Mli>  S  D&’1  J^l1  ' 


Proof 

For  £  =  1  the  result  follows  trivially  by  choosing  V  identically 

zero.  The  proof  proceeds  by  induction.  Let  7,  =  (T  }  be  a  connected 

Lh  i  i=l  ^ 

subset  of  £  and  let  (b  },  be  a  set  of  numbers,  with  £  b  =  0  ,  l  >  1  . 
“  i  i=l  i=l  ^ 

Select  T  €  ^  so  that  =  £^\{T)  is  connected  (it  is  easy  to  see  that 

this  is  always  possible);  to  simplify  notation  we  shall  assume  that  the 

n  i  h  h  h 

numbering  of  ^  is  such  that  T  =  and  that  T ^  and  T ^  ^  share  a  c 

edge.  We  define 

b^  ,  1  <  i  <  l  -  2 


ommon 


b. 

l 


b^,  ^  b^  ,  i  -  £  -  1  , 


and  use  the  induction  hypothesis  to  construct 


V  €  ?(4|>°<v">  .  pi‘1  •°<r) 


with 


(4.8a) 


V*V  =  0  at  all  vertices  of  , 
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(4.8b) 


V*Vdx=b^  ,  1  £i  <  t  -  1  ,  and 


l-l 


(4.8c) 


l£l!  i.Q(I")  s  DU-1)|‘'1  K!  • 


Let  W  be  the  field  constructed  in  Lemma  4.2  corresponding  to  the  triangles 
T ^  ^  and  ,  and  set 


(4.9) 


V  =  V  +  b^W  , 


where  V  and  W  are  interpreted  to  be  zero  outside  ^  and  ^(jT^ 

respectively.  This  V  clearly  satisfies  (4.7a)  and  (4.7b);  from  (4.9),  (4.8c) 
and  (4.6d)  it  follows  that 


(4.10) 


,  £-i  ,  , 

IY.iI  .  .<  D(£-l)h  1  l  | b  |  +  Dh  [  b  ^  | 

i,o(^)  i-i  1  l 


D£h_1  l  jb  |  , 


i=l 


(remember  that  D  at  all  points  in  this  lemma  is  the  same  constant  as  in  Lemma  4, 
This  completes  the  induction  argument.  | _ | 

Remark  4 . 2 

Based  on  (4.7c)  we  immediately  conclude  that 


(4.7c') 


hlli.flCj',  s  Dt3/2h_1  <  I  Ibjp) 

‘•h  i=l 


.2.1/2 


it  is  this  estimate  that  shall  be  used  later  on.  . _ , 


A  simple  rescaling  of  Lemma  2.5  in  [18]  leads  to  the  following. 


-  *- 
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Lemma  4 . 4 


Let  T  be  a  single  triangle  of  £  ,  and  let  <f>P  be  a  polynomial  of 


degree  <  p  such  that 


=  0  at  the  three  vertices  of  T  and 


dx  =  0  . 


There  exists  a  field  V 


P+1 


of  polynomials  of  degree  £  p  +  1  satisfying 


(4.11a)  VP+1  =  0  on  3Th 

(4.11b)  V-V^1  =  <J)P 

(*.iic)  l|vp+1||  1>Th  *  c(p+i)K||$p||  0>Th 

with  constants  C  and  K  that  are  independent 
of  ,  h  ,  p  and  $P  . 

In  this  lemma  we  have  again  used  the  fact  that  satisfies  a  minimal 

angle  condition. 

Lemmas  4.1  through  4.4  give  rise  to  a  local  construction  of  a  right 
inverse  for  the  divergence  operator.  We  give  the  details  of  this  construction 
with  particular  boundary  conditions;  this  result  shall  prove  useful  in  our 
proof  of  Theorem  5.1. 

Proposition  4 . 1 
Assume  that 


R(Ih)  *  6  >  0  , 


where  Rfl^)  is  the  measure  of  singularity  introduced  in 

independent  of  h  .  Let  J.  =  {T.}.  ,  denote  anv  subset 

‘-h  l  i=l 

be  any  element  of  P^p^’  )  ,  p  >,  3  that  vanishes 

vertices  of  .  Assume  that 


(3.9),  and  6  is 
of  ,  and  let 

at  all  boundary 


m 


24 


•P  dx  =  0  , 


for  any  connected  component  si"  of  T.(V')  .  There  exists  V  € 

-  h  —  4-  r-> 


)  X 


P[p+1]’°(i;)  such  chat 
(4.12a) 


V«V  =  p  in  sl(£')  ,  and 


(4.12b) 


Mi.ikD  s  c  pK  £3/2ll*!l0,fi(r >  • 


with  constants  C  and  K  that  are 
independent  of  ^  ,  h  ,  p  and  4>  • 

Proof: 

We  shall  without  loss  of  generality  restrict  our  attention  to  the  case  that 

H ([^)  has  only  one  connected  component.  Lemma  4.1  in  combination  with  Remark  4.1 

shows  how  to  construct  V  (  pt 3 ] , 0 ( J * )  x  p(3],0^c')  wit^ 

—1  u h  h 

<t>  -  V.V  =0  at  all  vertices  of  ^  . 


Lemma  4.3  applied  with 


bi  =  Lh  ^  4,~'7  *Y,1 )  dx  ,  l{i  i  l  , 
'i 


yields  €  p^l*0^)  *  P^’^(£h)  such  that 


-  7-(Vl+V2)  =0  at  all  vertices  of  ,  and 


(4>-V* (V^+V2) )dx  *  0  for  any  T  6  ^  . 


The  problem  is  now  completely  localized,  and  applying  Lemma  4.4  triangle  by 
triangle  we  find  V3  €  P^p+i^,0(^)  x  P^p+1^’°(^h)  ,  satisfying 


i 


3jj  "i 
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$  -  V-(V1+V2)  =  V»V_3 


i. e. ,  the  field 


V  =  v x  +  v2  +  V 3  C  P[p+1J,0(^)  x  P[p+1],0(^) 

has  the  desired  property  (4.12a).  It  follows  directly  from  Lemmas  4.1  and 
4.4  that 


(4.13) 


'-l''l,8(J') 


'o  Ml') 


Since 


i— 3II 1  n(TV  c  pK  (il*!i  v-  +  ^  Hi  suyV 
J  i,wvV  o.n(L)  j-i  J 


(^-V-V  )dx|  f  C  h(||*||  +  ||V  II  ), 

r  1  0 ,tJ  1  l.T" 


the  estimate  (4.7c)  shows 


Iloilo, n<D  5  c  t3/2(ll»Ho.fi(j')+ll*llli,n<r)) 


A  combination  of  (4.13)  and  (4.14)  yields  the  estimate  (4.12b)  for  V_  .  | _ | 

The  previous  argument,  with  minor  changes,  provides  proofs  of  both 
Proposition  3.1  and  Proposition  3.2.  Note,  however,  that  for  the  estimate 
(4.12b)  to  be  valid  for  $  €  ptpl>  ^(£^)  and  corresponding  V  €  p^P+^*®(£^)  x 
p[p+U,0(£h)  we  have  to  require  that  R(£h)  5  6  >  0  ,  independent  of 
and  h  (this  latter  is  the  reason  we  use  Proposition  4.1  and  not  the 
corresponding  version  of  Proposition  3.2  in  our  proof  of  Theorem  5.1).  If 
y,  is  taken  to  be  all  of  ,  then  t  'v  0  (h  ^)  ,  and  the  estimate  (4.12b)  reads 


IWl<0  5  c  •,Kh'3ll*llo,n 
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i.e. ,  the  local  construction  does  not  immediately  give  a  bound  for  a  right 
inverse  which  is  uniform  in  h  . 


i 


f 
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5 .  The  main  theorem 

As  announced  earlier  the  main  focus  of  this  paper  is  to  estimate  the 
norm  of  a  right  inverse  for  the  divergence  operator.  Our  estimate  is  the 
central  part  of  the  following  theorem. 

Theorem  5.1 

Let  V,  ,  0  <  h  ^  1  ,  be  a  quasiuniform  family  of  tr iangul at  ions  of 
the  polygonal  domain  *2  ,  and  let  p  be  an  integer  >  3  .  Assume  that 

R(J^)  is  >  0  ,  6  independent  of  h  , 

where  R(J^)  is  c^e  measure  of  singularity  introduced  in  (3.9).  Then 
V-(P[p+1]*°(Ih)  x  P[p+1]’°(Ih))  =  P[p],-1(Ih)  , 
and  there  exists  a  linear  operator 

:  P[p]’_1(Ih)  -  P[p+1]’°([h)  x  P[P+1]’°(j;h) 


such  that 

(5.1a) 

V-(Lh<P)  =  <P  v<f> 
P 

(5.1b) 

II  lN  II  *CpK 
p  i,n 

M0,« 

with  constants 

C  and  K  that  are  independent 

of  h  ,  p  and 

■t  • 

Note:  The  first  part  of  Theorem  5.1  is  simply  a  restatement  of  Proposition  3.1. 
Also  note  that  the  assumption  R(£^)  i  5  >  0  does  not  rule  out  the  presence  of 
singular  vertices,  it  merely  prevents  the  nonsingular  vertices  from  becoming 


too  close  to  singular. 
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Lj! 


Since 


V.(P[P+11*°([  )  X  p[p+l].°(y  ))  c  PIp],_1(L) 


it  is  well  known  that  the  statements  of  Theorem  5.1  are  equivalent  to  the  so- 


-1, 


called  inf-sup  condition  (c  =  C  ) 


(5.2) 


Plp,’'1(j>h) 
j,  V*V  <t>  dx 


sup 

V 


cp~K  Mlo.n 


1^1,1 


with  the  supremum  taken  over  V_  €  p^P+^>^(V  )  x  P^p+^^’^(V  )  (cf.  [2]). 


We 


-h'  '  v<-h' 

shall  make  use  of  this  fact  in  the  case  p  =  3  of  our  proof.  The  proof  of 
Theorem  5.1  relies  heavily  on  the  analysis  of  [IS],  but  an  added  new  element  is 
the  localization  procedure  which  has  certain  similarities  to  the  macro-element 
concept  found  in  [3,14];  however,  our  triangulations  are  quite  arbitrary, 
except  for  the  assumption  of  quasiuniformity. 


Lemma  5 . 1 


There  exists  a  constant  C  such  that  for  any  given  positive  integer  k  and 
h  sufficiently  small  (how  small  depends  on  k)  it  is  possible  to  partition  £ 


into  a  disjoint  union  of  connected  subsets  *  1  $  m  $  M(k,h)  with 


(5.3a) 


each  subset  containing  at  most  Ck  triangles, 


(5.3b) 


each  ft 


(m) 


of  radius  /kh 


interior (  containing  a  ball 

^h 
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Proof : 

Let  ,  1  $  m  <  K(k,h),  be  those  vertices  of  the  uniform  lattice, 

with  sidelength  2(/k+l)h,  that  lie  in  q  and  lie  at  least  a  distance  Jk.  h 
away  from  3ft  .  Let  D^m')  denote  the  open  disc  of  radius  //k  h  ,  centered  at 


4m)  •  All  triangles  of  that  intersect  will  be  assigned  to  the 

subset  *  thus  ensuring  that  (5.3b)  is  satisfied.  At  this  point  the  sets 

are  connected,  mutually  disjoint  and  each  contains  at  most  Ck  triangles. 

It  is  now  easv  to  distribute  the  remaining  triangles  of  £  among  the  ,  1° 

h  n 

such  a  way  chat  their  individual  connectivity  is  preserved,  and  they  still 
satisfy  (5.3a)  (possibly  with  a  larger  constant  C). 


(m) 


Remark  5.1 

Based  on  Lemma  5.1  we  may  immediately  conclude  that  for  h  sufficiently 
small  (how  small  depends  on  k5  it  is  possible  to  partition  ^  into  a 
disjoint  union  of  connected  subsets  ,  1  £  m  S  M(k,h)  satisfying 


(5.4a) 


each  subset 


y(ra) 

^h 


contains  at  most 


k 


triangles. 


(5.4b) 

The  constant  c  is 


each  ft^m^  =  interior (  ,  .Tilft) 

h  T^(m) 

independent  of  k  and  h  . 


contains  a  disc  of  radius 


cvic 


Lemma  5 . 2 
Let  k  be  a 
1  {  m  (  M(k,h)  be 
1  S  m  5  M<  k ,  h)  and 


positive  integer.  For  h  sufficiently  small,  let  » 

the  partition  of  £  introduced  in  Remark  5.1.  For  any 
any  constant  b  ,  one  can  find  €  P^’^(J^m^)  such  that 


b 


and 
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(m) 


0,0. 


(ra)  ' 


C(/kh) 


-li 


|U(m)||  ,,  5  C(/kh)-1(L+(v^h)~1)  |b| 


Proof : 


,(m) 


From  (5.4b)  we  know  that  there  exists  z  €  such  that 

—  n 


D  (c»£h)  c  ^m)  , 


where  D^r)  is  the  open  disc  of  radius  r  centered  at  £  .  Selecting  £  to  be 
the  origin  and  rescaling  by  c*^th  we  obtain 


D0U)  c  0 


(m) 


where  0^™^  is  the  translated,  rescaled  image  of  .  Let  'l''™*  be  the 

triangulation  of  m  corresponding  to  ;  it  is  then  possible  to  construct 


o  (m) 


•  (m) 


€  1 ,0(^-(m) )  satisfyj:ng 


(5.5a) 


=  0  on  the  boundary  of 


>  (in) 


(5.5b) 


.  (in) 


ip  dx  =  1  ,  and 


(5.5c) 


l*l,i.a«  5  C  • 


The  function 


(m) 


(x)  -  b  (c»kh) 


x-z 


c  *ich 
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°fll  0  (m) 

Is  an  element  of  *  (£,  )  that  satisfies  the  requirements  in  this 


Lemma. 


We  are  now  ready  for  the 


Proof  of  Theorem  5.1 


Consider  the  case  p  =  3  ;  we  shall  verify  that  for  h  sufficiently 
small  and  for  any  <J>  €  p^»  ^(7  )  there  exists  W  €  ptM.O^y  j  x  p[4]*0^r  ) 


(5.6a) 


(5.6b) 


s  1/2ll*H0.s  “d 


IHII  1>0  *  c  • 


It  follows  immediately  from  (5.6a-b)  that 


/  V  *V  <J>  dx  /  7  *W  <fi  dx 

SUp  -  :>  - 

v  mu  1>3  mu  1>a 


Jn<))2dx  -j^  ( <J>—V -W)  (J>  dx 


— '1,£2 


^  iisu1.;c,,*,,o-a  • 


1 . e. ,  the  inequality  in  (5.2)  holds  for  p  =  3  .  According  to  the  comments 
made  earlier  this  proves  the  theorem,  in  the  case  p  =  3  ,  for  h  sufficiently 
small.  For  p  =  3  and  large  h  the  theorem  follows  directly  from  the 
constructive  proof  of  Proposition  3.1,  discussed  at  the  end  of  section  4. 

The  construction  of  W  proceeds  in  several  steps. 


I 


Step  1:  Using  Lemma  4.1,  with  £  =  ^  and  p  =  3  one  finds 


Sl « pI3K0(s:h)  -  p,j,’udh) 


133,0. 


such  that 


(5.7a) 


V-V  =0  at  all  vertices  of  ^  , 


(5.7b) 


'!— ! II 1 , $  c 


Step  2:  Let  be  the  disjoint  partition  of  introduced  in 

Remark  5.1.  Let  $^ra')  €  be  the  function  constructed  in  Lemma  5.2 


corresponding  to 


(♦-V-V^dx  , 


and  define 


f>(x)  *  *(m)(x)  for  x  €  ,  1  k  m  k  M(k,h) 


It  follows  from  Letnma  5.2  and  (5.7a-b)  that 


(5.8a) 


p  -  7-V  -  |  =  0  at  all  vertices  on  t! 

boundaries  of  ,  1  $  m  £  M(k,h)  , 

h 


(5.8b) 


.  ,  (<j>-V-V,-<jj)dx  =  0  ,  1  k  ra  <  M(k,h)  , 

Am)  — i  — 


(5.8c) 


IIJIl0fQ s  c  Ikll  o  »  and  ll  *  'I  i.;-2  -  cfl+(,^h)  1)IUII0> 


Step  3:  The  function  t  -  7-V  -  1  i  P^’  ^  )  vanishes  at  all 

boundary  vertices  of  -.(^m^)  ,  1  i  m  i  M(k,h)  ;  it  also  satisfies 


J  (m) 
Uh 


( <f>-V  - V1-<f> )  dx  =  0  . 
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n  I  n  f  JJJ  j 

We  apply  Proposition  4.1  with  ^  =  to  this  function  for  each 

1  <  m  <  M(k,h)  .  By  composition  of  the  individual  solutions  we  get 


satisfying 

(5.9a) 


—2  € 


pi 4  ] 


<L) 


V  ’V_2  =  <P  -  V  ■  in  f!  ,  and 


(5-9b)  llv2ll1>0  ;  C  b3/2ll*ll(M  • 

Step  4:  Finally  we  shall  construct  a  field  ^  )  *  p[4],0^  ^ 

such  that 


(5.10a)  *  Oq^^*1)  IUIIq^  ,  and 

(5.10b)  112, II,.,  s  C  ll»ll0-fl  . 

In  combination  with  (5.7b)  and  (5.9a-b)  this  leads  to 


(5 . 11a)  lit  -7.W||0^  S  ^(h+vi:-1)  |U|lo,a  ,  and 

(5-ub)  ll-H  i.n  *  c  k3/2H  ♦  Ho.o  • 


where  W=  £v.  €  P141,0^)  x  P[4)’°([h)  .  If  k 
j=l  J 

large  and  h  is  sufficiently  small,  then  we  obtain 


is  chosen  sufficiently 


C0(h+/k-i)  <  1/2  , 

and  (5.11a-b)  therefore  verifies  the  existence  of  a  field  W  with  the 


properties  (5.6a-b). 
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1 


The  construction  of  is  based  on  an  approximation  argument.  Let  $ 

be  a  function  satisfying 


(5.12) 


(5.13a) 


(5.13b) 


A4>  =  $  in  ft  ,  with 


I*  «2.q  *  c  H»  «o.a  a”d 


*  II  3,(1  5  °  II*  HlJJ ' 


(note  that  we  do  not  specify  any  boundary  condition  on  3ft  ,  and  this  is 
what  makes  it  possible  to  obtain  (5.13a-b),  although  3ft  is  not  smooth). 


Let  V-  €  (T  )  x  pl^3»0(£  \  5e  an  approximate 


ion  to  V4>  in  the 


sense  that 


(5.14a) 


(5.14b) 


f  Ch  11*113,0  “ad 


Mi,o  5  c  ; 


(5.12)  and  the  estimates  (5.13b),  (5.14a)  then  lead  to 


•v3l|0i!J  -  IMMV«-V3)||0j( 


*  Ch  « 


s  Chi  c 


l,ft  * 


so  that  by  virtue  of  (5.8c) 


“♦-’•Mo.n  ss'w/ir-1)  IM0ifl 


The  remaining  inequality  (5.10b)  follows  immediately  from  (5.8c)  ,  (5.13a) 
and  (5.14b)  . 


EEL  J'-Jm 
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This  completes  the  proof  of  Theorem  5.1  in  the  case  p  =  3  . 

Let  p  be  an  arbitary  integer  5  4  .  Given  $  €  P^’  )  it  is 

h  r* 

possible  on  each  triangle  T  of  /,  to  find  a  quadratic  q  with 


(5.15) 


q  =  <p  at  the  three  vertices  of 


T*1 


dx  ,  and 


T11 


1  v.  II  e  $  C  h  supU(x)  |  <  C  pK  ||<fr  || 

r*  o.t*  x€T  o.r1 


for  K  >  , 

(in  (5.15)  we  have  used  the  Sobolev  Imbedding  Lemma  and  a  Bernstein-type 
inequality,  cf.  [18]).  Define 


q(x)  =  q  (x)  for  x  ►  f' 

~  v  - 

then  q  6  P[2,’_1(V  )  c  Pl3i’_1(V  )  ,  Sin,e 
i-h  —  *-h 

we  conclude  that 


(5.16) 


lkll0>n  f  c  pk  !!♦  I o , • 


T4' 


)  .  F rom  (5.1  )) 


and 


(5.17) 


o.ii  5  c  p  !lo,n 


Due  to  our  method  of  construction 


<)>  -  q  *  0  at  all  vertices  of  ][  ,  and 
|  =  ®  on  trlan8les  of  ^  . 


I  I 
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We  may  now  apply  Lemma  4.4  separately  on  each  triangle, and  by  piecing  together 
we  get 


vL  e  P[p+11,0(£h)  x  P[p+11’°(Ih)  , 


with 


(5.18a) 


V-V  =  4>  -  q  in  (2  ,  and 


(5.18b) 


|Vi«  l>a  £  C  p-  0>„ 


s  c  PK+K  ||*  || 


0,(2 


r  3  1  ___  ^ 

Since  q  €  Pl  ’  (£  )  we  may  use  this  theorem  in  the  case  p  =  3  (which  has 


already  been  verified)  to  find 


v2  €  P[4]’°(£h)  x  P[41’0(Th)  , 


such  that 


(5.19a) 


V-V  *  q  in  !2  and 


(5.19b) 


'-2 1*  1,(2  '  C  Hqll  0,(2 


$  C  p  ||4> 


'0,(2  ’ 


in  the  last  inequality  we  used  (5.16).  Defining 


V  -  vL  +  v2  , 


the  theorem  follows  directly  from  (5.18a-b)  and  (5.19a-b)  in  the  case 
p  5  4  .  This  concludes  our  proof.  | _ | 
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I  < 


)  » 


The  proof  presented  above  immediately  carries  over  to  the  case  of  homo¬ 
geneous  Dirichlet  boundary  conditions,  except  for  the  construction  of  4>  and 
V  .  We  need  an  additional  result  concerning  the  invertibility  of  the  diver¬ 
gence  operator  with  homogeneous  boundary  conditions.  The  following  lemma  is 
proven  in  [1];  the  method  of  proof  relies  heavily  on  the  characterization  of 
trace  spaces  for  function  spaces  on  polygonal  domains,  as  found  in  [8],  Sobolev 
spaces  H  (»J)  ,  0  f  s  ,  with  noninteger  indices  are  defined  by  interpolation; 

11*1!  _  0  denotes  the  norm  on  HS(il)  . 

Lemma  5 . 3 

Assume  that  all  internal  angles  at  corners  of  the  domain  Q  are  less 
than  2 tt  .  Suppose  that  <j>  €  Hs  (ft)  ,  for  some  0  <  s  <  1  ,  with 


$  dx  ■  0 

for  all  connected  components  of  .  Then  there  exists  U  €  HS 

such  that 

V  *IJ  =  $  in  ft  , 

U  =  0  on  9ft  , 


and 


^U.ft*  C  »*»lsffl  ’ 


*— '■  i.«  *  c  «♦»(>.<* 


with  C  independent  of  <t>  . 


38 


Let  $  be  as  introduced  in  step  2  of  the  proof  of  Theorem  5.1;  | 

1/2 

clearly  lies  in  H  (ft)  ,  and  it  has  integral  zero  on  each 

3/2  0 

connected  component  of  ft  .  Let  6  H  (ft)  fl  H(ft) 

be  the  field,  corresponding  to  $  ,  which  is  defined  by  Lemma  5.3.  If  ( 

°  f  1 1  o  °  r  1 1  o 

r  ’  x  '  ’  (£^)  is  an  approximation  to  in  the  sense  that 


IU3-V3I 


l,ft 


<  Ch 


1/2, 


-3"  3/2, ft 


and 


~~3  1 , ft  fC'l-3''l,ft 


then 


(5.20a) 


'V  -3^0, ft  =  HV'(V^3)|io,ft 


5  ‘IV^Ill.ft 


$  Ch1/2||U 


3 11  3/2 , ft 


<  Ch 


1/2  1 


1/2,  ft 


and 


(5.20b) 


Ml,ft  5  C' 


'0,ft 


Due  to  (5.20a),  (5.8c)  and  "logarithmic  convexity"  of  the  Sobolev  norms  it 
follows  that 

Hi-V.V3|l0>n  5  C(h+^r*1)1/2||4.||0  _  ; 
from  (5.20b)  and  (5.8c)  it  follows  that 
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1,0 


< 


For  k  sufficiently  large  and  h  sufficiently  small  therefore  has 

the  same  properties 


| II o( s  7 


and 


'-3' 


1.0 


$  C 


'0,ft 


as  the  field  constructed  in  step  4  of  the  previous  proof.  furthermore 

vanishes  on  3ft  and  hence  it  may  be  used  in  a  construction  of  a  field  with 
homogeneous  Dirichlet  boundary  conditions.  The  rest  of  the  proof  proceeds  as 
before,  thus  completing  our  verification  of 


Theorem  5.2 

Assume  that  all  internal  angles  at  corners  of  the  polygonal  domain  ft  are  less 
than  2tt  .  Let  £  ,  0  <  h  $  1  be  a  quasiuniform  family  of  triangulations 
of  ft  ,  and  let  p  be  an  integer  5  3  .  Assume  that 

O 

R(^)  ^  j  0  ,  6  independent  of  h  , 

O 

where  R(£^)  t*ie  measure  singularity  introduced  in  (3.10).  Then 

v.(pfp+1J’°(£h)  x  P[p+1],0([h))  =>  pfpJ.-1(^h>  f 
and  there  exists  a  linear  operator 

Lp  .  p[p]'-l(£h)  PlP+1J»°([h)  x  P[P+1J’°(^h) 


such  that 


(5.31a) 


V-Up4>)  -  ♦  ^  €  P[p],_1(Ih) 


40 


(5.3ib)  *  c  PK»*llo,n 

with  constants  C  and  K  that  are 
independent  of  h  ,  p  and  <J>  . 


Remark.  5 . 3 

Theorem  5.1  and  5.2  may  directly  be  used  to  show  that  minimization  of 
the  displacement  energy  of  two  dimensional  plane  strain  linear  elasticity  over 
the  space  of  continuous  piecewise  polynomials  of  degree  p  +  1  ,  p  >  3  , 
is  an  accurate  numerical  approach.  On  a  quasiuniform  family  of  triangulations 

O 

(with  R(Ih)  or  R(£h)  bounded  away  from  0)  it  leads  to  approximate  solutions 
that  coverge  at  optimal  rate  in  h  and  at  arbitrarily  close  to  optimal  rate 
in  p  ,  uniformly  with  respect  to  Poisson's  ratio  (cf.  [12,19]).  Theorem  5.1 
and  5.2  thus  disprove  a  conjecture  made  by  the  second  author  in 
Remark  3.2  of  [19];  it  was  conjectured,  based  on  numerical  evidence, 
that  the  h-convergence  rates  would  never  be  optimal,  uniformly  in  Poisson's 
ratio.  However,  the  numerical  experiments  referred  to  (cf.  [16])  were  all 
for  polynomials  of  degree  p  +  1  ,  p  <  3  ,  l.e. ,  exactly  the  case  the 
theorems  here  do  not  cover,  and  they  are  not  characteristic  of  the 
behaviour  for  p  ?  3  . 
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6.  A  basis  for  the  divergence  free  space 

In  many  applications,  it  is  of  interest  to  work  directly  with  the  null- 
space  of  the  divergence  operator  acting  on  ptP+2],0  *  p[p+l]»0  ^or 
,d[P+1],0  p[p+l],CL  A  ,  .  ,  , 

(r  <  r  )  .  As  observed  in  section  3  the  curl  operator  maps 

p[p+21,l(j  ^  (respectively  p^ P+- ^ ^ J ' ) )  onto  this  nullspace  (provided 
ft(£')  is  simply  connected).  Thus  a  basis  for  the  nullspace  can  be  obtained 
from  one  for  P^  (or  piP+2]Aj  ,  a  basis  for  pf P+2],l  wag  given  in 

[10].  We  shall  extend  slightly  that  work  here  to  construct  a  basis  for 
piP1^-  >  _  Our  method  of  proof  is  to  verify  the  dimension  formula 

(6.1)  dira(Plp+2],1(^'))  =  |p(p-5)T  +  (2p-l)EQ  +  3VQ  +  a  , 

°[d+21  1 

and  in  the  process  exhibit  this  many  linearly  independent  functions  in  PlF^  J> 
(these  functions  form  a  subset  of  the  basis  given  in  [10]);  T  here  denotes 
the  number  of  triangles  of  £  ,  Eq  ,  VQ  denotes  the  number  of  internal  edges  and 
internal  vertices  of  J  respectively  and  o  is  the  total  number  of  singular 
vertices  of  £  .  The  polygonal  domain  ft(£  )  is  assumed  to  be  simply  connected. 

The  operator  V*  maps  the  space 

P[P+l]’0(f)  x  ptP+1l»°(^') 

into 

p(p],-l(£')  . 

The  nullspace  of  V»  is  isomorphic  to 

and  it  thus  follows  that 
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(6.2) 


dim(Ptp+11’°  x  Ptp+1]’°)  _  dim(P^p+21,1) 
<  dim(Pfp},_1)  . 


The  first  and  the  last  of  the  dimensions  in  this  formula  have  already  been 
computed  to  be  (p-l)pT  +  2pEg  +  2V^  and  j(p+2)(p+l)T  -  a  -  1  respectively, 
i.e.,  based  on  (6.2)  we  get 


dim(P 


( P"b2 1  »  1  /  V  '  s  \  „  1 


(l  ))  5  jp(p-5)T  +  2pEQ  +  2V0  +  a  -  T  +  1  . 


Since  T  -  E  +  V  =  1  and  V  -  =  E  -  Eq  ,  this  implies 

(6.3)  dim(Plp+2l*1(J'))  5  (p-5)T  +  (2p-l)EQ  +  3VQ  + 


The  inequality  (6.3)  proves  half  of  the  identity  (6.1),  and  it  thus  remains 


to  verify  that 


(6.4) 


ditn(P 


[p+2],l,Y'<n  „  1 


(I  ))  *  p(p-5)T  +  (2p-l)E0  +  3VQ  +  a 


In  [10]  it  is  shown  that 


(6.5)  dim(P[p+2]’1(V'))  =  I(p+3)(p+4)T  _  (2p+5)EQ  +  3VQ  +  aQ 


through  the  construction  of  a  purely  local  basis  for  this  space.  Among  the 
corresponding  nodal  values  are 

(a)  the  value  and  and  derivatives  at  each  vertex, 

(b)  the  value  at  each  of  p  -  3  distinct  points  in  the 
interior  of  each  edge, 

(c)  the  (edge)  normal  derivative  at  each  of  p  -  2  distinct  points 
in  the  interior  of  each  edge. 


i 
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The  remaining  nodal  values  are  more  complicated  to  describe,  but  for  vertices 
on  the  boundary  of  fUj')  they  do  include 

(d)  one  cross  derivative  (i . e .  for  each  vertex  on  the  boundary, 

select  adjacent  edges  e^  and  meeting  there  and  take 

the  e^  ,  e^  cross  derivative  at  that  vertex) , 

(e)  the  second  edge  derivative  for  all  the  edges  meeting  there. 

For  functions  in  the  no<}ai  values  in  (a)-(c)  corresponding  to 

vertices  and  edges  on  the  boundary  of  £  must  vanish;  by  a  simple  count 
we  get  that 

(6.6)  3(V-VQ)  +  (2p-5)(E-E0) 

nodal  values  must  vanish.  The  second  derivatives  along  the  boundary  edges  at 
vertices  on  the  boundary  (e)  must  also  vanish,  and  give  rise  to  2  vanishing 
nodal  values  per  vertex.  Finally,  if  we  pick  e^  or  e^  in  (d)  to  be  one 
of  the  boundary  edges,  it  is  clear  that  this  produces  one  additional  nodal 
value  that  must  vanish  for  functions  in  P  p  ’  Q  )  .  In  combination  with 

(6.6)  we  get  a  total  of 

(6.7)  3(V-VQ)  +  (2p-5)(E-E0)  +  3(V-VQ) 

*  2p(E-EQ)  +  (V-VQ) 

vanishing  nodal  values.  Using  (6.5),  (6.7)  and  the  fact  that  E  +  Eq  *  3T 
and  E  -•  Eq  ■  V  -  Vq  we  thus  obtain 

(6.8)  dim(Plp+21,'1(r>)  *  p(p-5)T  +  (2p-l)E0  +  3VQ  +  o 

+  ((V-VQ)  -  (o-o0))  . 


44 


The  right  hand  side  of  (6.8)  is  exactly  as  desired  in  (6.4)  except  for  the 
additional  term  (V-V^)  -  (o-Oq)  ;  this  term  is  always  nonnegative  and  it 
equals  the  number  of  nonsingular  boundary  vertices .  In  order  to  verify  (6.4) 
it  therefore  suffices  to  find  one  nontrivial  linear  constraint,  for  the  nodal 
values  corresponding  to  each  nonsingular  boundary  vertex,  which  must  be 
satisfied  by  functions  in  p^P"1^  1  ^ (£ ’ )  .  a  constraint,  that  is,  which  is 
not  already  counted  in  (6.7). 

Let  Xq  be  a  boundary  vertex  and  let  the  triangles  ,  angles  0^ 
and  edges  e^  meeting  at  this  vertex  be  numbered  consecutively  as  shown  in 
Fig.  7. 


*o 


Fig.  7.  Boundary  vertex 

Let  denote  the  directional  derivative  in  the  direction  parallel  to  the 

edge  e  .  There  is  a  simple  relationship  among  all  the  cross  derivatives  of 
$  €  »*•(£')  '  namely, 
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(6.9)  sece.  3  3  ($  T  )(x~)  =  -  sec9.  .  3  3  (<p  T  )  (0 

l  e.  e . , ,  T|T.  —0  i-1  e.  ,  e.  ‘7.  .  —0 

i  i+1  l  l-l  l  i-1 

2 

+  (cote  +cot9.  . ) 3  ($|T  )(Xo)  , 

i  l-l  e .  ( .  ,  — u 

i  i-1 

2  $  i  f  r  (see  [10]  and  also  [5]).  Summation  of  (6.9)  with  alternating 
signs  yields 

r  i  2 

y  (-1)  (cote  +COte.  .  )  3  (4>|r  )0O  * 

,  i  i-1  e  'I .  .  — 0 

i=2  i  i-1 


sece  3  3  ( $  |  T  )  (2Sq)  +  ("D^ece  3  U  It  )  0O  . 

l  el  e2  /i  -u  r  er  er+1  -0 

For  *  €  P(p+2]’l(y:)  both  3  3  (4>IT  )  and  3  3  ($|T  )  must  vanish 

‘-h  e^  '7^  er  er+2.  r 

at  Xq  and  we  thus  arrive  at  the  constraint 
r  .  „ 

(6.io)  l  (-i)1(cote.+cote  )s^  (<j>!t  ) (x.)  *  o  . 

i=2  1  1-1  ei  i-1 

At  any  nonsingular  boundary  vertex,  r  is  at  least  2  and  cote^  +■  cote^_^  ^  0 
for  some  1  ,  so  that  (6.10)  represents  a  nontrivial  linear  constraint  among 
the  second  edge  derivatives,  which  is  not  counted  in  (6.7);  this  completes  the 
proof  of  the  identity  (6.1).  At  the  nonsingular  boundary  vertices  the  expression 

r  i  2 

l  (-1)  (cote  +cote  )3  ($|T  )(x~) 

i-2  i  i-1 

can  be  used  as  a  nodal  value  for  P^p+2^’1  in  place  of  one  of  the  second  edge  derivative* 
(one,  for  which  cot 9^  +  cot  9^  ^  +  0)  .  Using  these  nodal  variables  we 
obtain  a  basis  for  plP+2  ]*!■(£')  directly  from  the  basis  for 
by  deleting  members  corresponding  to  the  aforementioned 

2p(E-E0)  +  2 (V-VQ)  -  (o-oQ) 

vanishing  nodal  values. 
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Remark  6 . 1 

In  Che  case  U  is  not  simply  connected,  one  finds  that  the  nullspace  of 
V*  is  tlie  curl  of  the  subspace  in  p(p+-l'l  consisting  of  functions  that 
are  constant  on  each  component  of  312  ,  and  whose  normal  derivatives  vanish 
on  HI  .  This  space  lias  a  natural  basis,  and  its  dimension  exceeds  (6.1) 
exactly  by  the  number  of  components  of  ill  .  Using  the  corresponding  Euler's 
formula  we  can  thus  extend  our  combinatorial  proof  of  Proposition  3.2  to  domains 
that  are  not  simply  connected,  j — ^ 


I 
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